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ABSTRACT 



We report the discovery of the first variable extreme horizontal branch star in a globular cluster (oj Cen). The oscillation uncovered 
has a period of 114 s and an amplitude of 32 mmags. A comparison between horizontal branch models and observed optical colours 
indicates an effective temperature of 31,500+6,300 K for this star, placing it within the instability strip for rapidly oscillating B 
subdwarf s. The time scale and amplitude of the pulsation detected are also in line with what is expected for this type of variable, thus 
strengthening the case for the discovery of a new subdwarf B pulsator. 
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1, INTRODUCTION 

Subdwarf B (sdB) stars are hot (20,000 K < T^^ < 40,000 
K), compact (5.2 < logg < 6.2), evolved objects found on 
the extreme horizontal branch (E HB), a high t emperature ex- 
tension of the horizontal branch dHebei]|2008l) . Over the last 
decade, a small subset (~ 5%) of them have been found to ex- 
' hibit rapid luminosity variations with typical peri ods of 100-200 
, s and amplitudes on the mmag scale (see e.g. iFontaine et al.l 
■ l2006l for a recent review). Also known as EC 14026 stars af- 
ter the prototype (iKilkenny et al.|[T997|) . the rapid pulsators are 
found among the hotter sdB stars in a well-defined instability 
strip from 29,000-36,000 K. Pulsation driving is believed to 
occur through the action of a classical A:-mechanism associated 
with a local overabundance of iron peak elements in the driv- 
ing region, which in turn reli es on radiative levitation processes 
[ (ICharpinet et al.|[T996[ Il997|) . The modelling of these stars has 
' been very successful, to the point where their non-adiabatic pul- 
sation properties are understood not only qualitatively, but also 
quantitatively. Indeed, full asteroseismological analyses leading 
to a precise determination of the targets' fundamental parame- 
ters have now been car ried out for 12 out of at least 35 known 
EC 14 026 stars (see e.g. Ivan Grootel et a DI2QQ8I: ICharpinet et al.1 
1200 8' for recent results). 

It is generally accepted that sdB stars are the progeny of red 
giant branch stars that underwent significant mass loss around 
the time of the core helium flash, leaving them with very thin 
hydrogen-rich envelopes. These stellar structures do not expe- 
rience the asymptotic giant b ranch phase, instead evolving as 
so-called AGB-Manque stars (iGreggio & Renzini|[T99Q) . how- 
ever we still lack firm constraints on the evolutionary chan- 
nels) followed. A number of different scenarios involving e.g. a 
common envelope phase, stable and unstable Roche lobe over- 



flow, o r the merger of two H e white dwarfs have been pro- 
posed ('Han et al.| ]2002l l2QQ3l) . but remain to be tested obser- 
vationally. More recently, it has been suggested (iHanl l2QQ8h 
that the low fraction of cl ose binaries iden tified among clus- 
ter EHB stars (^ 4% from Moni Bidin et al | |20Q8.) compared to 
the field pop ulation (^ 60% ac cording to iMaxted et al.l [2QQll: 
^ 40% from iNapiwotzkil l2QQ6h can be explained naturally in 
terms of an age effect. In the case of an older population, the 
white dwarf merger scenario forming a single sdB star be- 
comes more important than common envelope evolution pro- 
ducing an sdB star in a close binary. It has also been suggested 
that cluster EHB stars might be the progeny of the "hot heli um- 
flasher" scenario ( Castellani & Castellani 1993; D'Cruz et al.l 
119961: ICasteUani et al.l l2QQ6l: iMiUer Bertolami et al.l |2008). In 
this scenario, red giant stars approaching the tip of the RGB 
lose a significant fraction of their envelope due to a vio- 
lent mass-loss event. The ensuing red giants have a higher 
mass than the limit for core helium ignition, and undergo a 
He-core flash either whe n approaching, or along the white 
dwarf cooling sequence dBrown et al.l l2QQll) . Another work- 
ing hypothesis is that cluster EHB sta rs are the progeny of 
the heliu m-enhanced sub-populations (D'Antona et alj 120021: 
iLee et al.l i2 QQ5 : D' Antona & Ventura 2007) recentlv d etected in 
a few globular clusters (iBedin et al.ll2QQ4l:lPiottoll20Q8l) . 



Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile (proposal ID 
081.D-0746). 



This is where asteroseismology becomes important: since 
the different evolutionary scenarios leave an imprint on the fun- 
damental parameters of the resulting sdB population, contrasting 
the predictions from population synthesis models with the cor- 
responding asteroseismic parameter distributions can be used to 
distinguish between the proposed format ion channels. First at - 
tempts in this direction have been made (iFontaine et al.ll2Q08h , 
however the number of EC 14026 stars analysed up to now is 
too small for quantitative conclusions. Moreover, the sample of 
known sdB pulsators has so far been confined to the Galactic 
field population. While several attempts have bee n made to de- 
tect variable EHB stars in globular clusters (e.g. iCatelan et al.l 
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Fig. 1. CMD of OJ Cen based on our SUSI2 w-band data com- 
bined with the WFI V-band data. The 52 selected subdwarf B 
star candidates are marked by large black points, and the posi- 
tion of the pulsating sdB star candidate is indicated. Error bars 
on the right display the mean colour and magnitude errors as a 
function of brightness. 



I2QQ8I: iKaluznv & ThompsonI l2008t iReed et al.l l2006l) . no con- 
vincing EC 14026 star candidates were uncovered. 

Rapidly pulsating subdwarf B stars in globular clusters are of 
high interest in the context of an asteroseismic characterisation 
of the EHB population in these systems. The first necessary step 
in this direction is obviously the detection of pulsating subdwarf 
B stars in globular clusters. This is why, when given a couple 
of free hours at the end of a night of observations dedicated to 
a diff'erent programme, we decided to monitor a field in oj Cen 
for rapid variability among the EHB population. We present the 
outcome of this mini- survey in what follows. 

2. OBSERVATIONS 

We obtained 2 hours of time-series photometry of a 5.5' x 5.5' 
field in the southeast quadrant of oj Cen with the SUperb Seeing 
Imager (SUSI2) mounted on the NTT at La Silla, Chile, on 14 
March 2008. A t/-band filter was used to minimize field crowd- 
ing, and 3x3 binning was chosen in order to reduce the overhead 
time to 16 s. Combined with an exposure time of 20 s, this re- 
sulted in a cycle time of 36 s, short enough to detect the rapid 
oscillations expected in EC 14026 type stars. 

After bias- and flatfield-correcting the images, photometry 
was performed using DAQ PHQT/ALLSTAR and ALLFRAME 
(IStetsonll 1 9871 1 1 99 ill 1 991 . Assuming a Moff'at analytical func- 
tion and using the task ALLSTAR, we first estimated an an- 
alytical point- spread function (PSF) for each frame by select- 
ing about 50 bright, isolated stars uniformly distributed over the 
chip. In order to obtain a global star catalogue where all im- 
ages reside in the same coordinate system, we used the task 
DAOMATCH/DAOMASTER ( Stetson 1994) to put together the 
192 images obtained for the field. As a reference catalogue, 
we employed UBVI-band photometry of oj Cen gathered with 



the Wide Fi eld Imager (WFI) on t he 2.2-m ESO/MPI telescope 
at La Silla (' Castellani et al.l 1200 7). which completely overlaps 
the SUSI2 field. Finally, we performed simultaneous PSF-fitting 
photometry on the SUSI2 data frames using ALLFRAME. The 
seeing on the individual images ranged from 0.76'' to 1.3", cor- 
responding to FWHM values between 2.9 and 5.4 pixels, while 
the photometric accuracy at u ~ 18.5 mag ranges from ~ 0.1 
to ~ 0.15 mag. The final star catalogue of the field monitored 
features -20,000 stars, of which 52 were selected as subdwarf B 
star candidates on the basis of brightness (16<w<18.5), colour 
(-2 < (u-V) < -0.8), photometric accuracy ([Su, ^Vlmean < 0.05 
mag), sharpness (s < 1), and separation index (sepj^^y > 10). In 
Fig. [1] we show the colour magnitude diagram (CMD) obtained 
by cross-correlating our SUSI2 w-band data with the correspond- 
ing V-band measurements from the WFI catalogue. Note that the 
data are not calibrated in magnitude; this is however of no rel- 
evance for the present study since we are interested only in the 
relative brightness variation of the targets. 



3. RESULTS 

Selecting only the 52 EHB stars indicated in Fig. [H we com- 
puted the light curves with respect to the mean w-band mag- 
nitude, incorporating corrections for airmass and seeing varia- 
tions by scaling the SUSI2 magnitudes to a WFI t/-band ref- 
erence magnitude. Fourier transforms were then calculated for 
each light curve in the 1-15 mHz range, appropriate for detect- 
ing the pulsations expected for the EC 14026 pulsators. We dis- 
play two selected pairs of light curves and the corresponding 
Fourier transforms in Fig.[2l The lefthand plot is particularly in- 
teresting since it features a strong peak well above the adopted 
Act detection threshold at a period of 114 s (8.75 mHz) and an 
amplitude of -30 mmag, which is compatible with EC 14026- 
type pulsations. The secondary peak is clearly an artefact of the 
data as it corresponds to the Nyquist frequency (or a period of 
72 s, twice the sampling time), and is encountered in the Fourier 



^ The separation index quantifies the degree of crowding 
dStetson et al. 2003). 
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Fig. 3. Finding chart for the variable EHB star discovered in oj 
Cen. 
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Fig. 2. Left Panel: Light curve (top) and Fourier transform (bottom) for the pulsating EHB star discovered in co Cen. The relative 
u magnitude indicated for the light curves is offset by around +0.25 mags with respect to the absolute magnitude from the WFI 
catalogue. The horizontal dashed line indicates the Act detection threshold. Right Panel: The same as in the left panel, but for a 
comparison star deemed to be constant. 



transform of most of the targets monitored. On the other hand, 
the stronger 114s luminosity variation is unique to the star dis- 
played, and it can be seen clearly in the phase folded light curve 
displayed in Fig.lH The righthand panel of Fig. [2] features a typ- 
ical EHB star from our sample that is judged constant on short 
time scales, with only the 72-s (13.89 mHz) artefact clearly vis- 
ible. 

The variable star uncovered has J2000.0 coordinates 
a= 13^27^ 11.77^. and d=-AT3T29.Q'' . A finding chart based on 
the combined SUSI2 image is shown in Fig. [3l It can be seen 
that the candidate pulsator is located close to two other stars 
of similar brightness, with a separation of ~ 1.5 times the av- 
erage F WHM. In fact, space -based ACS imaging of the same 
region (ICastellani et al.ll2QQ7h reveals a faint companion to our 
target that is not resolved by SUSI2 and is therefore included 
in the PSF computed for the target. However, the flux contri- 
bution of this companion is only ~1%, too small to potentially 
induce the >3% luminosity variation detected. A verification of 
the photometry shows no periodic change in brightness for the 
stars closest to our pulsator candidate, and moreover the resid- 
uals after PSF subtraction are negligible. We therefore believe 
that the 114-s luminosity variation observed is inherent to the 
star itself rather than constituting an observational artefact. 

Since the SUSI2 field overlaps with both the WFI and the 
ACS catalogues of oj Cen, we can determine the magnitudes 
and colours of the variable star in both filter systems. These 
are 5ac5 = 18.30, 7?ac5 = 18.55, if q^ac^ = 18.48, 16.81, 
5^^^=18.27, yw/=18.46, and 7^7=18.52. The location of the 
pulsator in the CMD obtained on the basis of both the WFI and 
ACS data is shown in Fig. [5l and reveals this star to be a typ- 
ical member of the cd Cen EHB population. We can therefore 
employ appropriate horizontal branch models in order to derive 
an estimate of the eff'ective temperature of our target. High S/N 
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Fig. 4. Phase-folded light curve for the variable EHB star, based 
on the 1 14-s periodicity detected. The sinusoid overplotted con- 
stitutes a least squares fit to the light curve and has an amplitude 
of 0.03 mags and a phase of 0.22. 



spectroscopy would yield more precise results, however we do 
not have such measurements at our disposal for the time being. 
In Fig. [5] we have overplotted the zero-age horizontal branch 
(ZAHB) and the exhaustion of central He-burning publicly 
available from the BASTI database (see, e.g. Pietrinferni et al.l 
12004 with an alpha enhancement [Qf/Fe]=0.4, helium abun- 
dance 7=0.248, and metallicity Z=0.002. We assumed a distance 
modulus yUo= 13.70+0.06 (D el Principe et a l. 2006) and redden- 
ing E{B - y)=0.1 1+0.02 ( Calamida et al.ll2005h . It can be seen 
that in both plots the pulsator is located close to current evolu- 
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Fig. 5. F625W, - F625W (left) and /, V - / (right) CMDs for selected stars of the a; Cen ACS and W FI catalogues. Details 

concerning the selection criteria for the star catalogues used in the plot are given in lCastellani et alJ (l2QQ7l) . The error bars on the 
right of each plot account for colour and magnitude photometric errors, and the number of stars in each catalogue is indicated. We 
overplot the ZAHB and the exhaustion of central He burning for the labelled composition. The candidate sdB pulsator is marked by 
an arrow and asterisk in both plots. 



tionary predictions within the theoretical and empirical errors. 
Given an estimated accuracy of ~ 10 % for the properties of 
the horizontal branch (HB) models and the slight dependency 
on the metallicity assumed (we averaged values for metallicities 
of Z=0.002 and Z=0.0006, appropriate for the metal-rich and 
metal-poor HB populations found in co Cen) we conservatively 
estimate Teff = 31,500 ± 6300 K for the variable star uncov- 
ered. This falls right into the 29,000-36,000 K instability strip 
for rapid sdB pulsators. The mass is estimated to be -0.5 M©, 
which is typical for a subdwarf B star. Given that the period and 
amplitude of the oscillation are also in line with what is expected 
for this type of pulsator, we believe that we have uncovered the 
first EC 14026 star in a globular cluster. 



for any sdB pulsator (the highest ampl itude pulsation is a round 
50 mmags for Balloon 09010001, see lOreiro et aP l2004l: typi- 
cal amplitudes are on the 1-5 mmag scale), and lower amplitude 
luminosity variations will be challenging to uncover given the 
faintness of the star and the crowded field. Nevertheless, we hope 
that the results presented here will constitute a first step towards 
the asteroseismology of globular cluster subdwarf B stars, and 
create the unique opportunity to provide quantitative constraints 
on their evolutionary history. 
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ular SUSI2 instrument scientist Alessandro Ederoclite for their support. Sadly, 
the observations reported here were among the last ever obtained with SUSI2 
since the instrument has since been decomissioned. 



4. CONCLUSION 

On the basis of 2 hours of SUSI2 time- series data, as well as the 
availability of WFI and ACS multi-colour photometry, we dis- 
covered what appears to be the first rapidly pulsating subdwarf 
B star in a globular cluster. Both the pulsation properties and the 
eff'ective temperature estimate are in line with what is observed 
for this type of variable in the field population; however, a longer 
time series of photometry and high S/N spectroscopy are needed 
to verify the results presented here. It will be particularly inter- 
esting to see whether more frequencies can be uncovered, since 
all known sdB pulsators show multi-periodic luminosity varia- 
tions. However, given the magnitude of the pulsator discovered 
and the shortness of the time series obtained, it is not surprising 
that only one periodicity was detected. Indeed, the period spec- 
tra of fast sdB pulsators are quite often dominated by one high- 
amplitude mode, and individual frequencies can be too closely 
spaced to be distinguishable in just one two-hour run. The am- 
plitude of the oscillation measured is among the highest found 
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